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A New Theory of Sources of Birkeland Currents 
by K. D. cole* 

Laboratory for Planetary Atmospheres 
NASA/Goddard Space Plight Center 
Greenbelt, Maryland 20771 

U.S.A. 


Abstract 

A new approach to collisionless plasma (Cole 1983) shows the existence of 
current orthogonal to ^ along the low latitude boundary layer of the magneto- 
sphere driven by electric field which is orthogonal to both ^ and the layer. 

In this case the relationship 

b 2 eEj_ 

p, + — - — -- = constant, 

8ti 8tt 

holds on a line orthogonal to B and the layer, where e is the dielectric 
constant of the plasma for electric fields orthogonal to Across the 
geomagnetic tail there flows a current in the direction of the dawn-dusk 
-electric field, and in this case a relationship 

B^ 

holds along a line orthogonal to _E and Bi^. Divergence of both these currents 
is shown to be a source of Birkeland currents. Also some of the boundary 
layer current is continuous with current across the tail. Electric currents 
of physically similar origin flow In interplanetary space, and when the 
magnetosphere interrupts them, additional Birkeland currents are driven. 


* On leave from La Trobe University, Bundoora, Victoria, Australia, 3083. 
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Infaroducfaltin 

In n rQCtsnfc, i^ipor a tjencsval fchoocy of oloGhfomai^notlo PiaXils ia plasmas 
in a qunsi-sliQatJy sb<5.to was proiinsovl (cola 1983). This bheory is now applicul 
bo tho problem of bUo soiu*oqs of olacbi'lc cuinrenb iresponsiblo fob goomagnobie 
tUsburbance. in blm flbsb pabb of bhe papeb, some iiloallsod configurabions of 
oloobbomagnoblG fioWs in plasmas aro oonoiilebeil. In bha soeomi panb, bhaso 
abo applied bo sbruGhvn*es in the goomagnetlG field and inbeiplanetaby spate ns 
pabb of a new bheoby of geomagnebio disbubhance. 

Thouglv many Gomponenbs of eai'lieb bhoobies of geomagnetic dlstubbance 
(Axfobd and Hines, 1961; Dungey, 1961; piddingbon, 1960,196?.; Cole, 1960, 

1961) remain valid today, none of them is onhirely satisfactory and all fail 
to be quanbibablvo on the precise mechanism of bransfer of energy from the 
solar wind into bhe magnetospbere and the abmosphere. Later, this author 
(Cole 1974) illustrated mocbanisms whex'oby solar wind plasma could enter the 
goomagneUic field. These included "grad B capture" and plasma flow across the 
magnebopanse caused by bangenhial aleobric fields and inertial drifts. 

The tangential electric fields could be caused either as the result of 
grad B capture on the day side of bhe magnobopause (Cole 1974) or as b)ie 
resulb of dynamo acbion in bhe Irlgh labibude ionosphere (cole 1976). 

In bhls paper an abbempb is made ho advance the theory of solar wind 
Intarachion with the magnetosphere invoking ne\«? understanding of currents 
generated by plasma flows (Cole, 1983) in interplanetary space, and in various 
regions of the magnetosphere and showing how the currents couple into the 
upper atmosphere. The units used are c.g.s. , except where otlierwise stated. 
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FiiUii conPlgmwbiona 

In bhliS i>apo^ fehn fehnocy of feho MaKwolL afei*oQt30S in n coiliaionimsiaj 
pUsuw (c^olo, I0{13) ia nppUo«\ ko aiacvnis fehn aouccoa oP olocbirlG cuirconb 
enapnnaibln Sot* nnomaqnnfelc ilinbucbanGQ. By ternnbing n Goillaionlissa plasma 
as a mo‘3inm wlUh spoctflahlo illninnbfic nnnsbanfe awl maqnofelc peemoahiliby ib 
haa boon ahovm bhab/ in n abaa<ly sbabo, 

(lUv p) E + (i K H) K n - ~ Ve “ "I i/vn » o, (i) 

whai'Q Vi « oloGbcic Pinltlf 

13 w waqnobio imlucbion, 

JB a inagnnbiG Pinlcl/ 
n « ~ mvignofeiq pncmoabiliby, 

wb.Bjre !i , denobe compcsnanbs pairallal ami pticpomllonlar I'oapecbivQly bo _b. 

, 2 
4im m e 

s 1 *»• I —• — r , 
n 

whobo Ug ^ luuubet* oP ions cm”^ oP spoGios b , 
and nig ^ moloGuIatr mass oP iona oP span ins s. 

Xb bas bonn shown nlsowhoro (Colo IBBS) bhab in plasmas such as Qjcist in 
bim. low labibntlQ baundacy layor atul inbotrplanebaty apace bha Povree denobed by 

1 O 

bhe boiau (div Jl ) E - ~ dqnnbion 1 delves slgnielcanb current. This 

Gurrenb is rePerred bo as dleleGtrlG cvtrrenb and arises principally because oP 
the large value oP e in, magnebised plasmas, in this paper divergence oP 
currents oP this kind in the geomagnoblc envlronmonb are invoked ns sources oP 
Birkeland currents observed in the magnobospbere. First, some simple models 
oP crossed olocbric and magnetic fields in plasnvns are presented, which define 
the new currents invoked. Then these models are applied to the plasma 
environment of the eiirth in apace ho estimate the strength of Birkeland 
cur rent a . 
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Prom equation (1) 

(div D) E “ 
- X 


X I n2 1 

2 Dx 2 ® 3x' 


8(B , ) 3(B , ) 

B ( z/u Jj/tL. 


3x 


3z 


3(b , ) 3(B . ) 

+ B.. . 0, 


3y 


3x 


(2) 


(div D) B - E ^ ^ - B ( - 

- y 2 1 3y 2 3y M x 3y 


3(3 > ) 3(B . ) 

2£/lL. _ - Y/y 


3x 


and 


(div D) E - 


3(B , ) 3(B , ) 

X B = 0 

^ ®z ^ 3z 8y ^ ' 


1 E 2 1 g2 3_ 1 _ 3 (ILjl/mI 

2 1 3z 2 3z p y 3z 


(3) 


°'Vu> 

3y 


3(B , ) 3(B , ) 

+ B.. 0. 


3x 


3 z 


(4) 


when e >> 1/ div d_ « eE^^ = — 1|^’ 

Consider now different simple cases of equation 1 which can be integrated and 
applied to interplanetary space or the geomagnetic field. 


Current of Type I ; 


Suppose B = (0, 0, B ) 
“ z 


(5) 


and E = (0, E (y), E ) 
- y •' z 

It is shown in Cole (1983) that 


( 6 ) 


P = 


1 


1 + 



This case (see Pig. 1), discussed earlier (Cole, 1983) shows that a current 
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flows in the k direction, of magnitude given by 


1 c 3 2 

'“y - ®"pj- 

z 


( 7 ) 


3 3 

Also if << Ey equation 3 intergrates, when ^ “ 0/ bo 


2 2 

Btip, + b - eE = constant, 
i 2 y 


(8) 


This corresponds to the circumstances of plasma flows and currents in the x 
direction. We investigate later the possibility that the flow in the low 
latitude boundary layer of the earth (Eastman, 1979jEastman and Hones, 1979) 
is like this* Also we expect that flows describable by this model exist in 
interplanetary space (see applications, later in the paper). 


Current of Type II ; 


In this case (see Fig. 2), 


suppose B = (0, 0, B ), 
— z 

and E= (0, E (x), E ). 
— y z 


with 


3y 


3x 


0, and « Ey equation 2 yields 


3 / z_x . „ 2 3e - 

^ ( — “) + E ^ • 

3x p y 3x 


( 9 ) 

( 10 ) 


( 11 ) 


Integrating equation 1 1 , yields 


2 2 3 2 

BttPj^ + B^ + eE^ - (E^ ) dx = constant. (12) 
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This caso d * appropriate to tilsousslons of flow of current across the 
"neutral" sheet in the geomagnetic tail under the influence of a "cross-tail" 
electric field. 

In this case the current in the y-directlon driven by electric field is 
given by 


. 1 c g. 2 j e 

Or B y 8x* 


(13) 


In tlie special case, when ^ E = 0/ equation 12 becomes 

2 2 

8irp, + D + eE = constant, 
i 2 y 

In the case of cold plasma p. =0/ and, 

^ 2 
2 C 

B ^ = constant;. 


(14) 


(15) 


where = C^/s = square of Alfven velocity 


(16) 


Equation (15) resembles the results obtained by Alfven (19S8) and Cowley 
( 1973) in discussing properties of neutral surfaces. While the general 

J 

solution of the problem of an infinite neutij^ sheet is contained in equation 
(14), the Alfven (1968) type solution is contained as a special case 
(corresponding to pj^ = 0), while the Harris (1961) type solution is another 
case, corresponding to By = 0, In the former case, a balance of forces is 
maintained by-V(eE )/8it and c ' j x B, In the case E„“ 0, a balance of j x B 
and “Vpj^ could be maintained. 

Case II is, of course, relevant to discussion of magnetic "merging," 
Equation 14 can be interpreted as implying the existence of static structures 
of opposed magnetic fields in which plasma is driven towards the merging 
region by an applied electric field Ey. The plasma is turned away from the 
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neubfal lino boCoiTG ceachlrig it genorabing Gui'cenh in the lUceeUlon oC By- 
The i'jabtlbioning oP on«rgy batwoon Uhormal energy of the plasma^ aa maniPesbefJ 
by |3 , magnetic energy ami elecbtlc field onecgy density/ is dependonb on the 
plasma density, or equivalently, the AlPven speed (see equation 14). clearly 
bhera is a limitless variety of possible partitions. Experiments looking for 
merging regions in the goomagnohic hail or ah the magnetopause sliould take 
into account all three terms of equation 14. It appears that a steady state 
neutral sheet can exist in which no merging takas place. Additional factors 
to those considered hero would appear to ho necessary In order to generate 
merging. 

A Variant of Case I: 


Assume n « 0, b„) 


(17) 


and ^ “ (0, E (y), 0) 


(13) 


D S 

Equation 3 then yields, assuming ~ « q t* •— 

u 5C u S 


(19) 


. X 1 ~ 1 n> 2 M 4, 1 L.f1t » n 


rd _wsdt*. 4 _ X / 1 

r~(eE ) E “"E r-t-B ^("-) - B 

Lgy y ■' y 2 y 9y 2 9y M x 3y 


Therefore, 




S , V, 2, 9 ,B , 

( SE ) V' (' ) ** 0 , 

By ' y 9y u ' 


and it follows by integration along the y-direction, that, 


(20) 


( 21 ) 


2 2 

bE ~ B “ 8irp, « constant. 

y i 


( 22 ) 
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This is essentially the same result as case I which shows that the components 
of B perpendicular to By act independently but similarly in relation to By. 


General two dimensional plasma; 


In this case, 


assume JB » (0,0, ~B^) 


(23) 


£= (Ejj, Ey, 0) 


(24) 


and ^ » 0 
02 


(25) 


Equation 2 new yields 




ay'-y' 


1 _ 2 8e 
9x 


1 8 , 1 . _ . z/\i _ 

+ 2 ® ^ ‘ 3x ’ ■ 


(26) 


ft a 1 2 9 E 

" 2 ® 3y 


From (27) 


1 29 , 1 , 

2 "x ^ 


- B 




[“(EE )1 E + (eE ^) 

’•3y y X 9x' X 


^(b!) 

9xS ’ 


z 9y 


1 2 9e 

2 ^y 9x 


= 0 . 


= 0 . 


(27) 


(28) 


From (28) 


r3 1 9 , 2, 9 ,B , 1 „ 2 9s . 

[k'"®x>l \ 37'"®y > - ■ 2 ®x 37 ' “■ 


(29) 


8 


♦ •'i 


since 7 X E » 0, equation 29 yields 
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|;(^) 


0, (30) 


and equation 30 yields 

b-\‘< 




(31) 


It follows that 


3k 




2 Lx 


1 eE 2 ^ B_] 

2 y y J 


2 «2 


3y 


1 r 2 B 

2 1 ww * ^ OXif ** 

L y 2 X y 


[eE^ 


] (32) 


No simple interpretation can be given at this time in this case of a general 
two dimensional plasma ami it would appear to need computer solutions to 
produce further erudition. The symmetry of equation 32 is, however, noted. 


Intermei^iiiate 7e and 7(1/y) 

Suppose a unidirectional electric field exists orthogonal to ^ and let 
there be gradients of e and y orthogonal to B, so that. 


= (0, 0, B^), 


(33) 


E == (0, Ey, 0), 


(34) 




(35) 


0 ^ 0 


(36) 


Equation (2) then yields, for the x~component , 


1 r. 2 3e . 1 _2 3 

2 ®y 3 


3(B 




( 37 ) 
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and for the y-component, 


r^~(eE ) “ E ~ 


B r^{B , ) ^ 0i 
"z9y'"z/p' 


(38) 


Equation 38 integrates to 


eE^^ - n^-8itpj^ w f(x). 


(39) 


while equation 37 yields 

eE ^ + B^+8np, = <{)(y). (40) 

y 

Once again it appears that computer calculations could be required to 
make physical interpretation more accessible in this case. 


APPLICATIOr^S 

We are now in a position to discuss the various current regions of the 

magnetosphere. First we consider the regions near the magnetopause. 

Traditionally/ currents at the magnetopause have been considered to be 

produced solely by a gradient of pressure pointing outwards from the 

magnetosphere. A Lorentz force jj. x B was considered to oppose the force due 

to this gradient of pressure (see Willis 1971 for a review). These currents 

are the chapman-Perraro currents. It is considered here that the traditional 

approach must now be modified. It will be seen to be Oi; limited application 

in portions of the magnetosphere only, viz./ over a small range of longitudes 

either side of the "nose" of the magnetosphere. In this region a negative 

pressure gradient acts from the solar wind into the magnetosphere. In vast 

regions of the magnetosphere/ and the low- latitude boundary layer associated 

with it/ there is a flow of solar wind plasma inside and roughly parallel 

to the magnetopause in which the speed decreases inwards (see e.g. / Fig. 3 

2 

taken from Eastman (1979))/ and in which the term eEj^ is a considerable 
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2 

fraction of p . In making this claim, we assume that eE^ is estimated by V 2 
(c,f., Cole, 1983). 

Boundacy Layer 

It is claimed by the present author that in the boundary layer the 
current system is different from what might be expected on chapman~Ferraro 
theory and the consequences are significantly different. One long recognized 
difficulty of the chapman-Ferraro approach is that it predicts no energy or 
plasma flow into the magnetosphere at all. Nor does it account for the 
neutral sheet in the geomagnetic tall. The chapman-Ferraro theory however 
appears applicable in a limited manner to some parts of the magnetopause as 
will be demonstrated. Let us now apply equations 7 and 8 of Case I to the low 
latitude boundary layer. 

Tlie Boundary Layer Current 

Consider an idealized low latitude boundary layer in which the flow of 
solar wind plasma is laminar and parallel to the boundary. Let the velocity, 
pressure and density be functions of distance from the magnetosheath. The 
magnetopause is considered to be a tangential discontinuity of bhe magnetic 
field. 

Figures 4a and 4b illustrate the situation in cross section in the dawn-side 
and dusk-side boundary layer. This configuration conforms to Case I. We 
realize from equation 7 that the dielectric current is sunward in the dusk 
side and anti-sunward in the dawn side. The presence of these currents has 
not hitherto been taken into account in the problem of solar wind interaction 
with the magnetosphere. Figure 5 illustrates the dielectric current in 
equatorial cross section. The current due to -Vp is not illustrated. 
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The dielectric current flows in the opposite direction to the classical 
Chapman-Ferraro currents which are caused by the gradient of plasma pressure. 
The dielectric current is a dielectric current of type I. The complete 
current orthogonal to ^ in such flows leads to equation 8 and is specified by 
equation 7. 

This theory calls into question the widespread view that the geomagnetic 
field is separated from the interplanetary medium solely by a current layer 
caused by the gradient of pressure of solar wind normal to the current 
layer. However the theory appears to be consistent with observations of the 
outer regions of the geomagnetic field where it interacts with the solar 
wind. Taking the ion plasma parameters for a typical dawn-dusk meridian in a 
sequence of steps through the plasma boundary layer (Eastman 1978) Table 1 is 
constructed 

TABLE I 

Plasma Magnetosheath PBL PBL PBL PBL Magnetosphere 

parameter 1234 

n(cm“^) 18 14 10 6 2 0.5 

V(km s"’') 300 280 200 120 50 <30 

i (kev) 0.5 1 2.5 4 6.5 7 

V2nmV^ 9.7 5.3 2.0 0.4 0.02 <0.003 

(kev cm“^) 

2 — 

p (= - nE) 6 9.3 16 16 8 2.4 

(kev cm"^) 

In this table E is the mean energy per particle. Therefore there is a 

substantial dielectric current to deal with. There is a decrease 
1 2 

of — p V (where Vg^ is the plasma velocity and p the density) from the 
solar wind into the boundary layer. 
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Ife Is evidonfe from Tablo 1^(a) fchnb the Gomponenfc off ourrenb in fcho Xayor 

auo Uo bUo grfuUonb of (gE^) is Of con\pv\b«hle ntagnibudo bo bhose clue bo 

grndlonbs of presauro (b) bho compononb duo bo gradienb prossuire odd a to bhab 

duo bo bbo oubQir aocblon of blio plosma boundacy Xayob (o) in bbo 

innoL* aocblon of bho loyor bbo GVP'ironb duo bo gundionb of prosauro is in bbe 

dti'QGbion of bho cXassionX Cbnpnuvn-Forrnro cuuconb. 

Tbo sbroncjbb of bbo dloXocbbic cui'ronb in bbo boundary Xayor can bo 

— 1 

osbimabod from ociuabion 15 assuming E “> c 

Then, por em nXong B, inbogrntod blurougli bbo Xayor ib ia osbimabad by 

(cm ^ ~ p V ^ B ^ . (<n) 

Tbo bobnX dioXoGbeio curronb (J.^) ia given by inbograting above and boXow bho 
oquaborial piano whorovor bbo boundary layer oxisba. Tbia ia nob known, bub 
wibbin a facbor of 2 ib ia osbimabed bo bo 10 Rj^ above and beXovf bbe 
Gquabox'ial piano* 


So 



'\j 


O t ^ * s • 


( 42 ) 


lUabXo TI sbowa values of for a varieby of values of B and n^ 

where p “ n„ m„ and m,. ia assumed « mass of n abom. 

S 9 8 
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V 

^sw 

n 

(amps) 

cm“^ ) 

( kms~ ^ ) 

(y) 


10 

400 

40 

4.S X 10^ 

10 

400 

10 

1.8 X 10*^ 

5 

800 

40 

9 X 10® 


These currents are geophysically slgniPicanb. They ©Kerb an inwards force on 
the magnetosphere and would conbrihube bo a decrease of the geomagnetic field 
earbh“ward of them, and an increase of magnetic induction outside the layer in 
the solar wind. Together with the tail current (see later in the paper) those 
currents cause a perpetual "ring" current which should manifest itself at the 
earth's surface as a depression of the geomagnetic field. All currents in the 
oarth environment need to bo taken into account in the search for such an 
effect, including field-aligned currents, otherwise known as Birkelnnd 
currents. Bquahion 0 is derived on the assumption that the fields and medium 
parameters do not vary in the x direction. In the low latitude boundary layer 
we find that they do vary. This would mean that current generated locally 
according to eqviation 42 would have divergence, it is now suggested that tl\e 
current in the low latitude boundary layer finds continuity partly with 
Birkeland currents into the polar ionosphere in the am sector and out in 
the . pm sector. Later the continuity with electric current across the 
geomagnetic tail will be discussed. Consider the quantity 

div Li 10 div{^), (43) 
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we note bhat bhe boundary layer .inside the geomagnetic field allows solar wind 
plasma to expand along the geomagnetic field towards the ionosphere. It 
follows that (pB S decreases with distance from the "nose" of, the 
magnetosphere. Approximately we can say that 



M 


T 


T 


M. 


m T 


(44) 


where M^, Ap, B = mass content, volume, length, cross sectional area and 
magnetic field (at equatorial plane) in a tube of magnetic field. Given 
and ApB (tube flux) constant during drift along the boundary layer, div jXp is 
negative. It follows that Birkeland current must flow along the geomagnetic 
field to the pre-noon sector of the polar ionosphere to which field lines from 
the magnetopause are connected (see Fig. 6), there to be continued by iono- 
spheric current. By symmetry Birkeland current would flow away from the polar 
ionosphere in the post-noon sector, it is suggested here that the part of 
this dielectric current which remains after flow along the flanks of the 
magnetosphere becomes continuous with current across the tail (see later). 
Birkeland current must flow along the geomagnetic field to the pre-noon sector 
of the polar ionosphere to which field lines from the magnetopause are 
connected (see Fig. 6), there to be continued by ionospheric current. By 
symmetry Birkeland current would flow away from the polar ionosphere in the 
post-noon sector. It is suggested here that the part of this dielectric 
current which remains after flow along the flanks of the magnetosphere becomes 
continuous with current across the tail (see later in the paper). 

The strength of the current per cm of longitude into the auroral 
ionosphere can be estimated by 

j II (cm ^ ) = f div (45) 
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where f is the separation at the boundary layer of the lines of force which 
are 1 cm apart in longitude near the cusp at the ionosphere. Approximately f 
= 15. Estimating the length scale of variation of p/B around the boundary 
using equation 45/ a value of 10 Rg is adopted. Then, at the ionosphere, 


. . - 1 , 

3, 1 i ■ 


(46) 


where the parameter p and B are of course boundary layer values. 

With B = 40y, = 4 X 10^ cm sec"^, p = 10 x 1.6 x 10”^'* gm cm"’^, 

-2 _i 

j|l = 1.0 X 10 amp cm . Such currents are in the range actually observed 
over the polar ionosphere (Sugiura and potemra, 1976). 

Let us return for a moment to the Chapman-Ferrero currents* In the 

region around the noon-midnight plane of the magnetopause the gradient of 

pressure of the solar wind may be considered as the dominant non-magnetic 

force and reflection of solar wind particles would cause a surface current 

which is towards dusk, equatorwards of the geomagnetic cusps and towards dawn, 

polewards of the cusps (see Fig. 7). Outside of this region the term 
2 

eEy /8ti appears to play a progressively more significant role increasing up 
to a substantial fraction of p. This is suggested by the work of Eastman 
(1979), who shows that in the numerous low latitude boundary layer crossings 
reported generally is between one half and one third of p. In the 

present context we estimate 


12 2 
~ p V « eE /Stt. 
2 sw y 


(47) 


It is clear that the magnitude of Chapraan-Ferraro currents near the noon 
plane will dominate the dielectric current elsewhere. This is so because the 
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reflection of solar wind particles in that region is equivalent to a gradient 

of pressure which, integrated across the region gives l\n equivalent 

2 

pressure (p ) of approximately 2 pV . The region in which reflection is the 

sw 

dominant interaction of the Solar wind with the magnetosphere is here called 
the "median” strip (illustrated in Fig. 7). 

Grad p current exists completely around the magnetosphere in the layer 
with density given by 


B X Vp 
i (Vp) » — 


the boundary layer it is estimated by 





sw 


Integrated throughout 


where p^^ and are values in the magnetosheath rather than the boundary 
layer. In the magnetosheath, along the flanks of the magnetosphere, the flow 
is expected to be approximately laminar with no variation of and 

therefore no divergence of j^(7p). it is suggested that near the "nose", 
spatial variation of Pgy,/Bg^ may be such that the ratio Increases away from 
the nose. This would form field aligned (Birkeland) currents towards the 
ionosphere in the a.m. sector and opposite in the p.m. sector. Also the 
Birkeland currents would tend to cause a depression of the magnetopause 
surface in the "median" strip between the cusps (see Fig. 9). However, 
polewards of the cusps, the Birkeland current from the Chapman-Ferraro current 
would reinforce that from the dielectric current, producing a dawn-to-dusk 
electric field at ionospheric levels and at intermediate heights in the 
magnetosphere . 
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HistorioalXy, the increase of geomagnetic fiel<3 at the earth's surface at 
the commencement of some magnetic storms has been interpretecl in terms of an 
increase of the Chapman-Perraro currents at the magnetopause particularly on 
the day side. This explanation would appear to be qualitatively valid in the 
present theory however the effect is complicated by the diminishing effect of 
dielectric currents and Vp currents in the opposite direction to the 
conventional Vp current of Chapman and Ferraro. 

In the region of the cusp of the magnetosphere and the median strip 
polewards of the cusp it is conceivable that a mixture of Chapman-Perraro 
current and dielectric current could flow if there is penetration of solar 
wind plasma there with a component of velocity orthogonal to B_. in this "zone 
of confusion" one could expect considerable variability in these currents and 
associated Birkelaid currents. This is suggested as the origin of the 
irregular Birkeland currents observed at ionospheric altitudes at the site of 
the dayside cusp. (Ledley and Farthing, 1974). 


The "Distant" Cross-Tail Current 

Par downstream from the eartli at distances greater than about 15 the 
geotail is known to contain an almost neutral sheet and attempts to explain it 
have, been in two classes. The first, exemplified by the work of Harris (1963) 
employs a set of "hot" particles only and essentially produces the result, on 
a line orthogonal to the sheet. 


Birp 



constant. 


(48) 


The second approach exemplified by Alfven (1968), works with cold plasma, and 
essential], y produces the result 
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2 2 

B + eE * constant (49) 

z y 

In a sense both approaches are valid but the most general approach combines 
the two, and is embodied in equation 14/ viz.. 


2 2 
Birp + B + eE 

2 y 


constant. 


(50) 


where Ey is now the cross-tail electric field, B^ bhe geotail 'tiagnetic field, 
and p the pressure of the entrapped plasma perpendicular to B, . 

It is proposed here that the cross-tail electric field is created under 
"ordinary" conditions by the dielectric current in the flanks of the low 
latitude boundary layer of the magnetosphere. Any mismatch of the current 
down the flanks and current across the tail would be taken up by Birkeland 
current at their junctions discharging through the polar ionosphere. 

The current density across the tail from the dawn side to the dusk side 
is given by 


j 


y 


c E 2 ^ 

B 8x “ 8tt B y 3x 


(51) 


The component of j^. 


per unit of tail length, may be estimated as. 


jy (cm S 


5B 

2tt 


amps cm 


-1 


(52) 


So that, with B = 10 Y/ jy (cm“^) = 8 x 10“^ amp cm”^. For a tail of scale 
length 40 Rj^, the total cross tail current from this cause would be 2 x 10^ 
amps. 
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9y a similar argument made in the case of the divergence of boundary 
layer currents to obtain Birkeland current/ a case is now made for such 
currents derived from the divergence ofcross-tail current using equation 
(52). This equation gives the locally produced current in the tail. It is 
assumed that changes of local parameters across the tail in the dawn-dusk (y) ‘ 

ii 

I 

direction lead to divergence of which leads to Birkeland current. Firstly/ . ij 

ii 

let us note the symmetries of the terms in j . Both-8p/3x and-9e/3x are I 

^ ii 

' j! 

expected to point away from the neutral sheet/ and since B reverses sign I 

ii 

there/ divergence of jy will not change sign with change of position above or || 

below the sheet. However the divergence of jy would be expected to change f 

sign depending on whether the position is on the dawn side or the dusk side of 
the center line of the sheet/ if there is symmetry of the local parameters 
about this line. It is expected on average that at any given distance X from 
the neutral sheet/ |^/ r— and B would each increase in magnitude from the 

o X o X 

central meridian plane of the tail to the magnetopause. The divergence of jy 
would then, be such as to produce Birkeland current towards the polar 
ionospheres from the "dawn" half of the tail, and from the ionospheres in the 
"dusk" half. However great departures from this simple pictures may be 
expected because of the dependence of div jy on 3p/3x, 3e/3x and B. This 
perhaps accounts for the complexity of Birkeland currents near the "Harang" 
discontinuity. 

The divergence of tail current may be estimated using equation 52, by 

3y - « j,l (amps cm"''), (53) 

where W, L = half width, length of tail » “ separation in the 

tail of two lines of force 1 cm apart in longitude near the midnight auroral 
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belt. Therefore, for B » lOy and " 20, the Birkeland current per cm to 

each polar ionosphere would be « 1.2 x 10 amps cm“ \ Birkeland currents of 
this magnitude are observed in the night time auroral belt (potemra et al. 
1979). 

This distribution of Birkeland electric current would also be compatible 
with Sugiura's (1975) inference from observations of magnetic fields regarding 
the connection of Birkeland currents near the polar ionosphere with those from 
the tail. These currents would cause the magnetic field lines to be splayed 
outwards away from the earth between them and the neutral sheet and to be 
oppositely splayed between the currents and the tail mantles of the magneto- 
sphere (Fig. 8) . There would be a tendency for the field aligned currents in 
the tail to produce a component of B, orthogonal to the neutral sheet, between 
the neutral sheet and the mantles, pushing the magnetic field out at the 
middle of the tail mantles. This should cause a "fin" to exist in the 
magnetic field down the middle meridian of the mantles, (see Fig. 9). 

The electric field across the tail from dawn-side to dusk-side, will 
cause the drift of plasma towards the earth and the consequent energization of 
trapped plasma; production of partial ring current and Birkeland currents at 
places of divergence of the partial ring current (Harel et al. 1981). It is 
not the purpose of this paper to examine these internal processes of the 
magnetosphere but to discuss principally the external origin of the Birkeland 
currents, and the dawn-dusk electric field. 

Additional Currents from the Interplanetary Medium 

It is clear that equation 1 predicts significant electric currents in the 
interplanetary medium (Cole, 1983), and that such currents may find conduction 
paths through the geomagnetic plasmas if the magnetosphere encounters tliem. 
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The strength of these currents will now be estimateil an«l their interaction 
with the geomagnetic field discussed. The currents of interest are those 
associated with shears in the flov^ of interplanetary plasma and in neutral 
sheets. These are examples of types I and II treated in section 1. Both of 
these cases may apply to neutral sheets, in type I the electric field is 
perpendic'nlar to both B and the neutral sheet. In type ii the electric field 
is perpendicular to B but parallel to the neutral sheet. 

2 

In interplanetary space it is generally expected that p << so the 
current can be estimated in magnitude in both cases by 


j « 


c 


J 

Sir B L 


(54) 


where L is the scale of variation of a critical parameter e or E of the 

p 

system. When L > c 'E m./Be , the plasma drift velocity is cE /B (Cole 1983)/ 

X y 

so we can say, in that case. 


1/2nm,V c 
i s w 

B~L ' 


(55) 


where V_,. is the component of the solar wind orthogonal to B. Fig. 10 shows a 
range of values of j for a variety of values of Vg^, B and L. Assuming n = 10 
cm . There may be occasions, too, when currents due to pressure are 
significant. 

It is assumed that an interplanetary current, if Intersected by the 
magnetosphere, will find continuity across the magnetopause with field aligned 
current to the polar ionosphere and out again, or may cross the tail of the 
earth via the neutral sheet. This would be done by a build up of space charge 
at the boundaries to adjust the electric fields to give current continuity 
(see Fig. 11). This seems a reasonable proposition. 
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Suppose bhe egfecbive area of bhe magnebosphere accessible bo bha 
inberplanebary currenb is (15 if L > 15 R^. and 15 L R^ if L < 15 Rg, 

Then we see bhab for bhe sbrucbures characberized in Fig. 10/ may produce 
currenbs bhrough bhe magnebosphere of order 10 - 10 amp. Now observabions 

of magnebic bays (Cole 1962a) show bhab their durabions range from about 20 
mins, to about 200 mine, though some are outside these values. By assuming 
these bays are caused by changes to the magnebosphere/ionosphere system 
brought about by interplanetary currents flowing through the system/ it is 
inferred that structures of scale length greater than about 5 x 10^® cm have 
sufficient current to account for geomagnetic bays. 

Let us now consider the morphology of large scale interplanetary 
structures corresponding to types i and ii. Take a system of cartesian 
coordinates in which x is radially outwards from the sun/ y is azimuthal and z 
perpendicular to the plane of the ecliptic. Let the solar wind velocity be 
given by 


V 

— sw 





(56) 


Let 


For L » m f — / 


B = (B , B / B ). 
X y' z 


(57) 


E = (V B - V B , V B - V B / V B - V B ) , 
— yz zy zx xz xy yx 


O, -V B / V B / 
X z X y 


(58) 


and/ 


2 2 

BxB=V(B +B)/VBB/VBB. 
— — X X y xyx xzx 


(59) 
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The magnitude of the current j is then approximated by 


„ 2 . 2 . „ 2 , 

e V (B + B ) 

4 « sw y z 

Sir L Be ' 


C (B ^ + B 
y z 

Bit L B 


(60) 


where M is the ratio of the solar wind speed to the Alfven speed. 

In type i the direction of the current is given by ± JE x ^ depending on 
2 

whether V(eEj^ ) is in the same or opposite direction to E^(Cole 1976). In 

this case, ^ is in the •!' z direction if ^ is away from the sun and will be 

2 

in the plane of the ecliptic in the direction ± E x B if 7(eEj^') is in the ± z 
direction (respectively). 

In type II the direction of the current is in the direction of E for a 
2 

gradient of in the plane of the ecliptic and 2 will be in the ± z 

direction depending on whether B_ is away from or towards the sun. In this 
case we should expect current to flow through the magnetosphere from N to s or 
vice versa depending on whether ^ is towards or away from the sun, producing 
significant hemispherical asymmetries in magnetic disturbanc<^ as observed at 
the earth's surface. 

Generally for type I currents (i.e., corresponding to Fig. I geometry, if 
2 

V(eEj^ ) is in the direction of ± E , ^ is in the direction of ± JB x JB 

2 

respectively, and if L » mE c /Be, we can say that the direction of j is ± 
that given by equation 59. Whereas, for type II currents, if V(eEjJ“) is in 
tV\e direction of ± ^ x jB, is in the direction of ± E^ respectively, i.e., as 
speexfied by ± that given by equation (58). 

The amount of current converted to Birkeland current from the tail would 
depend upon the angle (a) between the current direction in interplanetary 
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space and the direction defined by the vector cross product of the geomagnetic 
axis direction and the normal to the geomagnetic neutral sheet. However 
Birkeland current from the day sectors of the magnetopause would be controlled 
more by the angle (6) between the current direction and the geomagnetic 
axis. The precise description of the amount of Birkeland current generated by 
divergence of interplanetary currents depends on the configuration of the the 
draping of interplanetary magnetic induction lines over the magnetopause. The 
detailed configurations of the plasma in this coupling process will be the 
subject of future work. 

Recently lijiraa and Poterma (1582) correlated a variety (20 in number) of 

interplanetary quantities with morning and afternoon Birkeland current 

densities and found that some are better correlated than others. One of the 

2 2 1/ 

best correlations is with the quantity [nV^^ (B^ + B^ ) s in 0/2 ]'2 where 0 

is the angle between the geomagnetic axis and the z-direction. This 
interplanetary quantity is clearly closely related ''mathematically'' to the 
current j defined by equation 60# although its origin in physical discussion 
is different. 

The theory of the present paper would explain a permanent Birkeland 
current system associated with divergence of dielectric current in the 
boundary layer and the geotail plus a variable component associated with 
interplanetary parameters via equation 60 and the angles a and 6. To test 
this in the manner of lijima and Potemra (1982) it would be preferable to use 
total Birkeland current per cm of auroral oval rather than current density. 

On the Question of Magnetic Merging 


The physics of case I and its variant (equation 22) are similar. 
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in caao I fchavo is no moi*gintj of magrvofclc field 

,U’, 

bQcauao the volociby of plasma is pfu*allol bo^^Useontlnaity. 

Type IX could dosocibe a ''tranqonbinl" dloGontinulfcy bub we see In this 
simple goomebuy bhnb bhofe is no merging, Additional factors, peirhapa time 
variation or eui*vatvu'e of JR or dissipation may be needed to be incorporated in 
the model to manifest merging. It is clear that oven without mervging 
geophysically significant electric current flow in the boundary layer and 
geomagnetic tail. 

The Bow Shock 

That part of the bow shock which is called a quasi-perpendicular shock 

resembles the conditions for type IX in the foregoing; while in the regions in 

which it is quasi-parallel, the conditions are those like type I, since the 

1 ? 2 

bhorronl pressure of solar wind plasma is ^ P "^sw" ^ clear that 

the electric cvxrrents in the bow shock are dielectric currants of the kind 
discussed in this paper, viz., typo II currents in quasi-perpendicular shocks 
and type I currents in qunsi-pax*allel shocks. The size of the currents in the 
bow shock are estimated to bo of the order of 10“ - 10 ' amps and represent 
another potential source of Birkeland currents from them to and through the 
magnetosphere to the ionosphere. This could liappen if there were field lines 
connected to the shock, at places of divergence of the dielectric current 
draping over the magnetosphere* 

Evaluation of Ring Current and Ionospheric Dynamo Currents 

Traditionally the ring current has been ascribed mostly to trapped 
energized magnetosphoric plasma on close goomagi'ietic field lines and 
manifested itself by a decrease of the geomagnetic field at the earth's 
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surface. The Chapman-Ferraco currenbs had the opposite eFCect, Given the 
validity of the new currents in space around the earth a re-evaluation of the 
component of depression of the geomagnetic field at the earth's surface due to 
conventional magnetospheric ring current will have to be done. Likewise a re- 
evaluation of the compression due to Chapman-Ferraro currents will need to be 
done, fhis would also affect the amount of geomagnetic variation attributable 
to dynamo action in the ionosphere. 

Ionospheric Dynamo Action 

The fact that neutral winds at ionospheric levels can produce significant 
electric fields and currents (Nagata et al. (1950); Fukushima, (1953); Cole, 
(1960, 1961) should be taken into account at all times. 

These fields and currents are particularly significant (1) when currents 
and fields applied from the boundary layer and tail are weak i.e., under 
conditions of relative geomagnetic "quiet", (2) as a non-negligible component 
of disturbances, because externally applied electric fields dj'ive ionospheric 
electric currents, which in turn accelerate the neutral gas (Cole, 1962b) by 
the Lorentz force (^ x E[) , otherwise known in the ionosphere as "ion drag" . 
This, in turn, produces an ionospheric dynamo electric field. (3) W^en solar 
wind fields are "turned off" after a period of activity the newly established 
winds in the ionosphere drive currents orthogonal to in the ionosphere which 
are in the opposite direction to those driven by the solar wind fields (Cole, 
1966). They will tend to set up electrostatic fields which cause convection 
of the magnetosphere also in the reverse sense to that caused by the solar 
wind fields. Moreover, the mismatch of electrostatic potentials, tended to be 
established this way, between the northern and southern hemisphere will cause 
field-aligned magnetospheric electric currents (cf. Cole, 1971). 
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Dlscusalon 

Having applied the unabridged Maxwell equations bo the discussion of the 

electrodynamics of plasmas (Cole 1983) new "dielectric" currents have been 

predicted in the magnetosphere and interplanetary space to be of geophysical 

significance. The low latitude boundary layer is the site of one such current 

in the direction ^ x ^ (type I) caused by an electric field (E ) orthogonal to 

2 

^ accompanied by a gradient of eE^^ in the direction of _E. Divergence of this 

current produces Birkeland current ho the polar ionosphere and the remainder 

is continuous with current across the geomagnetic tail. In the tail, the 

2 

cross-tail electric field is accompanied by a gradient of eEj^ perpendicular 
to ^ and ^ producing a current in the direction of E . In addition there is 
cross tail current caused by -V p . Divergence of the cross-tail field 
appears to be responsible for the Birkeland currents to the auroral region of 
the ionosphere which were inferred to exist by Sugiura (1975). 

In addition to these currents, dielectric currents of origin in 
interplanetary plasma structures may be conducted across the magnetopause 
through the magnetosphere by Birkeland currents, and also across the tail by 
alteration of the cross tail electric field. The paper does not address the 
problem of how interplanetary currents negotiate the bow shock bo enter the 
magnetosphere . 

This theory gives a new approach to the investigation of geomagnetic 
disturbance. It is of interest to note that the type II current was in fact 
invoked from a different standpoint by Alfven (1968). It is now seen to be 
one of a wider class of electric current hitherto overlooked in plasma 
physics. The currents come about when a plasma is subject to electric fields 
(charge distributions) and currents of external origin. The bulk of plasma 
physics appears to be concerned with completely internal processes and it is 
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usual to assume electrical neutrality exists except in special circumstances, 
thus ignoring the kinds of new effects discussed here. 

The theory so far developed using this new approach to collisionless 
plasmas (Cole, 1983) has been applied only to situations in which ^ is 
unidirectional. To develop the theory further one needs to express the "y" of 
the plasma in terms appropriate to general magnetic induction JB, i.e. to allow 
for effects of "parallel" pressure of plasma in addition to "perpendicular" 
pressure to which this paper is presently restricted. 
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Figure Captions 

Fig. 1 Geometry for type I current. 

Fig. 2 Geometry for type II current. 

Fig. 3 Example of parameters measured in low latitude boundary layer/ 
taken from Eastman (1979). The boundary layer is identified between the 
vertical full and dashed lines. 

Fig. 4 Schematic representation of velocity/ electric field energy 
density and dielectric current (type 1) in low latitude boundary layer. 
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Fig, 5 Scheme oE dielectric currents in low latitude boundary layer/ and 
tail. Currents due to Vp not illustrated. 

Fig. 6 ScJiematic representation oE Birkeland currents produced by 
divergence of currents in low latitude boundary layer and geotail. 

Fig. 7 Illustraton of the "median"' strip enclosed by dashed curve on day 
side magnetopause. 

Fig. 8 Illustration of "spalying" of geomagnetic field in tail due to 
Birkeland currents otherwise known as field-aligned currents (F.A.C.). 

Fig. 9 Illustration of dayside "valley" and night side "fin" on 
magnetopause surface. 

Fig. 10 Curves of measure of dielectric current in interplanetary 
structures assuming n m^^ = 1.67 x 10 g cm . L = scale length across B_ of 
structure. 

Fig. 11 Illustration of interplanetary currents impacting on the 
magnetosphere. They are an additional source of Birkeland currents inside the 
magnetosphere. 
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